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ABSTRACT 

The  existence  of  organized  structure*  in  tur- 
bulent shear  flow  has  been  the  subject  of  recent 
observational  discoveries  in  both  the  laboratory  and 
in  the  atmosphere  and  ocean.  The  recent  work  on 
modeling  such  structures  in  a temporally  developing, 
horizontally  homogeneous  turbulent  free  shear  layer 
(J.  T.  C.  Liu  and  L.  Merkine  1976  Proc.  Roy.  Soc. 
London,  Ser.  A,  Vol.  352,  pp.  213-247)  has  been 
extended  to  the  spatially  developing  mixing  layer, 
there  being  no  available  rational  transformation 
between  the  two  nonlinear  problems.  The  basis  for 
the  consideration  is  the  kinetic  energy  development 
of  the  mean  flow,  large-scale  structure  and  fine- 
grained turbulence  with  a conditional  average, 
supplementing  the  usual  time  average,  to  separate 
the  nonrandom  from  the  random  part  of  the  fluctu- 
ations. The  Integrated  form  of  the  energy  equations 
and  the  accompanying  shape  assumptions,  is  used  to 
derive  "amplitude"  equations  for  the  mean  flow, 
characterized  by  the  shear  layer  thickness,  the 
nonrandom  and  random  components  of  flow  which  are 
characterized  by  their  respective  energy  densities. 
In  general,  the  large-scale  structure  augments  the 
spreading  of  the  shear  layer  and  enhances  the  fine- 
grained turbulence  by  taking  energy  from  the  mean 
flow  and  transferring  it  to  the  turbulence  as  it 
amplifies  and  subsequently  decays.  The  maximal 
amplitude  of  the  large-scale  structure  is  attained 
by  the  Initially  most  amplified  mode,  however,  the 
relative  enhancement  of  the  fine-grained  turbulence 
is  achieved  by  both  the  magnitude  of  the  large- 
scale  structure  and  its  streamwise  lifetime.  Thus 
a greater  enhancement  of  the  turbulence  is  achiev- 
able by  the  lower  frequency  modes  which  have  longer 
straamwlse  lifetimes.  The  large-scale  structure 
can  also  be  controlled  by  increasing  the  initial 
level  of  turbulence,  which  would  render  its  decay 
more  rapidly. 


P = P + p 


NOMENCLATURE 


P.  P.  P' 


dimensionless  velocities  of  the  mean, 
large-scale  structure,  fine-grained 
turbulence 


dimensionless  pressure  of  the  mean, 
large-scale  structure,  fine-grain*: 
turbulence 
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dimensionless  coordinates  (x,  horizon- 
tal; y,  vertical) 

dimensionless  time 

average  which  produces  a mean  quantity, 
(the  time  average  in  the  spatial  prob- 
lem, horizontal  average  in  the  temporal 
problem) 

conditional  average 
mean  turbulent  stresses 

conditionally  averaged  turbulent 

stresses 


rij 5 <uiY  * Vj 


modulated  turbulent  stresses 


Re  Reynolds  number 

Other  quantities  are  defined  as  they  appear. 


INTRODUCTION 

Recent  reviews  of  laboratory  observations  have 
pointed  towards  the  omnipresence  of  large-scale 
organized  structures  in  turbulent  shear  flows,  with 
each  type  of  shear  flow  having  its  own  peculiar 
organized  structure  (1),  (2),  (3),  (4),  (5).  Sev- 
eral international  conferences  have  also  been  held 
in  which  such  organized  structures  in  turbulent 
flows  were  discussed  (6),  (7).  It  is  apparent  that 
turbulence  in  shear  flows  are  not  as  random  as  orig- 
inally thought  to  be  and,  in  his  Drvden  Research 
Lecture,  Roshko  (4)  reminds  us  that  Drvden  (6)  had 
pointed  cut  that  it  is  necessary  to  separate  the 
random  from  the  nonrandom  processes  in  turbulent 
flow.  This  view  was  also  emphasized  by  Liepmar.n  (9). 
Observations  strongly  suggest  that  the  large-scale 
organized  structures  are  important  to  the  mixing 
and  development  of  the  shear  flow  and  it  is  begin- 
ning to  be  recognized  that  the  study  of  their  cvna- 
mics  is  fundamental  to  anv  studv  of  turbulent  shear 
flows.  Such  studies  being  technologically  impor- 
tant to  combustion  and  reactive  flows  (7),  (101  and 
jet  noise  (11),  (12)  to  mention  a few.  Such  orga- 
nized structures  are  also  observed  to  occur  ir. 
regions  of  strong  local  shear  provided  bv  frontal 
systems  and  internal  waves  of  a much  larger  scale 


in  the  atmosphere  and  in  the  ocean  and  are  now 
thought  to  play  an  Important  role  in  the  spectral 
cascading  process  and  in  mixing  (13),  (14),  The 
observational  unmasking  of  the  large-scale  orga- 
nized structures  ir.  turbulent  flows  has,  at  least 
in  the  authors'  view,  given  turbulence  research  new 
and  exciting  lmpecus. 

We  shall  focus  our  attention  to  free  turbulent 
flovs  in  this  paper  in  which  the  mean  flow  velocity 
profiles  are  inflexional.  The  presence  of  the 
large-scale  structure  is  more  or  leas  attributable 
to  the  dynamical  instability  inherent  in  such  flows 
(9),  (12),  (15),  (16).  Although  the  large-scale 
organized  structure  may  be  deterministic,  It  is 
conceivable  that  it  is  subjected  to  random  forcing. 
Take  the  mixing  region,  for  Instance;  the  large- 
scale  structure  as  a consequence  of  upstream 
forcing  at  a given  frequency  will  propagate  down- 
stream at  a phase  velocity  characteristic  of  chat 
frequency.  Suppose  now  a slower  large-scale  struc- 
ture Is  first  set  up  and  subsequently  a faster  one 
Is  generated,  a "collision"  of  the  latter  with  the 
first  will  take  place  downstream  as  the  faster  one 
catches  up  with  the  slower  one.  It  is  not  entirely 
inconceivable  that  this  provides  one  of  the  expla- 
nations of  the  "pairing”  or  "vortex  coalescence" 
process  that  has  been  observed  (17).  This  essen- 
tially emphasizes  Che  importance  of  etudylng  the 
large-scale  organized  structure  under  well- 
controlled  rather  than  under  natural  conditions 
(9),  much  In  the  same  spirit  of  Scbubauer  and 
Skramscad  (18),  for  purposes  of  physically  Isolat- 
ing and  understanding  the  mechanisms  Involved. 

While  it  Is  generally  Important  to  understand 
the  coalescence  of  large-scale  organized  structures 
in  a turbulent  shear  flow.  It  seems  logical  that 
theoretical  understanding  of  the  Interactions  be- 
tween a single,  monochomatlc  large-scale  structure 
with  the  fine-grained  turbulence  be  obtained  with- 
out complicating  the  picture  with  the  more  compli- 
cated dynamical  effects.  On  the  other  hand, 

"vortex  coalescence"  Is  known  to  occur  even  with- 
out the  presence  of  fine-grained  turbulence  and 
contributions  towards  Its  theoretical  explanation 
can  be  made  from  purely  dynamical  considerations 
(19).  We  shall  consider  here  the  first  major 
problem  discussed  above,  that  of  the  Interaction 
between  the  large-scale  structure  and  the  dlspar- 
ately  fine-grained  turbulence  in  a free  turbulent 
shear  flow,  in  particular,  the  mixing  region. 

Our  aim  is  towards  physical  understanding  of 
the  lntersctlons  between  the  disparate  scales  of 
motion  which  lead  to  the  growth  and  development  of 
the  shear  layer,  rather  than  emphasizing  numerical 
prediction  or  eccuracy.  This  necessitates  approx- 
imate considerations  in  order  to  echleve  our  aim, 
one  which  may  be  helpful  towards  the  eventual 
incorporation  of  the  large-scale  structure  dynamics 
in  predictive  schemes  or  cowards  well-controlled 
experimental  unmasking  of  those  mechanisms  which 
contribute  to  the  understanding  of  the  development 
of  organized  structures  In  turbulent  shear  flows. 

To  proceed,  we  essentially  follov  earlier 
suggestions  (8),  (9),  (20)  and  theoretically  sort 
out  the  monochromatic  large-scale  organized 
structure  from  the  random  fine-grained  turbulence 
bv  a conditional  averaging  procedure  which  has  been 
described  and  realized  in  the  laboratory  elsewhere 


(21),  (22),  (23).  This  forms  the  starting  point 
for  the  subsequent  approximate  considerations  for 
the  study  of  simultaneous  interactions  between  Che 
mean  flow,  the  large-scale  structure  and  the 
fine-grained  turbulence.  We  have  reported  on  our 
previous  study  of  an  idealized  mixing  region  (24), 
in  which  the  mean  flow  consisted  of  horizontally 
homogeneous,  oppositely  directed  streams.  The 
shear  layer  and  Interactions  develop  in  time 
rather  than  In  the  straamvlse  distance.  Since  all 
mean  quantities  there  are  horizontally  homogeneous, 
there  is  no  streamwise  "diffusion”  in  the  mean  and 
this  comprises  one  of  the  natural  simplifications 
of  the  problem  (there  being  no  backwards  diffusion 
in  time  because  of  causality,  of  course).  Another 
simplification  in  Ref.  (24)  was  that  the  only  non- 
equilibrium  behavior  of  the  fine-grained  turbulence 
Is  due  to  Interactions  with  the  large-scale  struc- 
ture. The  large-scale  structure  there  develops  in 
time  and  Is  horizontally  periodic.  The  problem 
considered  there,  then,  was  the  temporal  problem 
and  there  does  not  exist  a physically  rational  one- 
to-one  transformation  to  the  streamwise  developing 
ehear  layer  or  spatial  problem  because  of  the 
causality  condition,  among  others. 

In  this  psper  we  shall  present  new  results 
pertaining  to  the  spatial  problem.  The  locally 
equilibrium  (with  the  mean)  behavior  of  the  fine- 
grained turbulence  Is  also  relaxed  so  that  in  our 
Integrated  energy  considerations  Its  evolution  is 
due  to  the  balance  or  inbalance  between  production 
from  both  the  mean  and  large-ecale  structure  sod 
viscous  dlssipstlon  rather  than  attributing  Its 
nonequilibrium  behavior  to  the  production  mecha- 
nism from  large-scale  structure  elone.  In  order 
to  fix  our  ideas  as  to  the  physical  picture  which 
might  emerge  from  our  analysis,  It  would  be  helpful 
to  recapitulate  the  Interesting  results  reported  by 
Binder  and  Favre-Harlnet  (23).  A well-controlled 
large-scale  structure  was  Imparted  upon  the  ehear 
layer  of  a round  jet  end  measurements  were  taken 
on  the  jet  axis.  In  the  absence  of  the  large- 
scale  structure  the  turbulence  level  on  the  jet 
axis  gradually  Increased  downstream.  With  the 
presence  of  the  large-scale  organized  structure, 
the  fine-grained  turbulence  level  Increased  at  a 
faster  rate  and  became  more  vigorous  with  the 
downstream  distance.  Simultaneously,  the  large- 
scale  structure  amplified  and  subsequently  de- 
cayed, as  if  having  extracted  energy  from  Che 
mean  flow  and  subsequently  transferring  It  to 
the  fine-grained  turbulence.  The  spreading  of  the 
jet  was  also  enhanced  by  the  presence  of  the  large- 
scale  structure.  Although  geometrical  effects  in 
the  experiments  certainly  play  an  important  role, 
the  qualitative  physical  situation  should  hold  for 
the  plane  mixing  region  that  we  are  to  consider 
here. 

FORMULATION 

The  expllclty  study  of  the  large-scale  organized 
structure  coexisting  with  random  fine-grained  turbu 
lence  necessitates  the  splitting  of  a flow  quantity 
into  three  components  consisting  of  the  mean  motion 
and  the  nonrandom  and  random  parts  of  the  fluctu- 
ations. The  usual  time  average,  over  a period  at 
least  chat  of  the  organized  structure,  produces  the 
equations  of  mean  motion.  A second,  conditional 
average  Chen  enables  one  to  sort  out  the  nonrandom 
part  of  the  fluctuations  from  the  random  part,  in 
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contrast  to  and  different  froa:  filtering.  Such  a 
splitting  procedure,  suitable  for  the  present 
study  is,  of  course,  not  new  and  has  been  discussed 
and  utilized  in  the  literature  (2),  (8),  (9),  (15), 
(16),  (21),  (22),  (23),  (25).  However,  the  first 
systematic  derivation  of  the  conservation  equations 
for  the  three  component  splitting  procedure  is 
attributable  to  K.  C.  Reynolds,  as  much  as  the  two 
component  splitting  procedure  was  to  0,  Reynolds. 

Ue  shall  consider  an  Incompressible  fluid. 

The  velocity  components  and  coordinates  are  made 
dimensionless  by  the  corresponding  characteristic 
scales  appropriate  to  the  problem.  For  instance, 
for  the  mixing  region,  the  characteristic  velocity 
is  the  free  stream  velocity  and  the  characteristic 
length  the  initial  boundary  layer  thickneaa.  The 
preaaure  is  made  dlmenslonleas  by  the  free  stream 
dynamic  preaaure  and  the  time  by  that  formed  from 
the  characteristic  velocity  and  length.  The  con- 
ditional average  la  deno_ted  by  <(  )>  and  the 
usual  time  average  by  ( ).  The  mean  flow  quanti- 
ties axe  denoted  by  capital  symbola,  the  large- 
scale  atructure  identified  by  the  supererlpt 
tilde  ' and  the  random  fine-grained  turbulence 
by  the  supererlpt  prime  ’ , The  conditional  and 
time  averagea  of  the  turbulence  quantities  are 
zero  by  definition.  The  conditional  average  of  a 
large-acale  organized  quantity  reproduces  Itself 
and  its  time  average  is  zero.  Because  the  two 
components  of  the  fluctuations  are  not  correlated, 
the  time  average  of  their  producta  is  zero.  The 
conditional  average  of  the  products  of  two  fine- 
grained turbulence  quantities  produces  a mean 
Reynolds  stress  and  a modulated  stress,  which  are 
crucial  to  the  interactions  between  the  random 
fine-grained  turbulence  with  the  mean  flow  and  the 
large-scale  structure,  respectively.  The  modulated 
turbulent  stresses  oscillste  at  the  frequency  of 
the  large-scale  organized  structure. 

Just  as  in  the  study  of  the  two-component 
0.  Reynolds  turbulence  problem,  for  all  practical 
purposes,  one  is  not  Interested  in  the  details  of 
the  turbulence  but  in  the  stresses  it  sets  up  and 
the  interaction  with  the  mean  flow  in  an  averaged 
sense,  in  the  explicit  considerstlon  of  the  large- 
scale  structure  one  is  also  not  interested  in  the 
details  of  the  fine-grained  turbulence  but  only  in 
its  modulated  stresses  and  its  mean  stresses  as  it 
interacts  with  the  large-scale  and  mean  motions. 
Finally,  the  large-scale  structure  interacts  with 
the  mean  flow  through  its  own  mean  stresses.  These 
can  be  best  illustrated  by  writing  down  the  momentum 
equations  for  the  mean  flow  and  the  large-scale 
structure 
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where  D/Dt  is  the  substantial  derivative  following 
the  mean  flow,  Re  is  the  Reynolds  number.  In  the 
absence  of  the  large-scale  structure,  or  if  it  were 
not  sorted  out  from  the  random  part,  the  problem 
reduces  to  the  usual  two-component  turbulent  flow 
problem.  In  the  absence  of  the  fine-grained  turbu- 
lence, the  problem  then  provides  the  starting  point 
of  nonlinear  hydrodynamic  theory  in  a laminar  flow 
(26).  It  is  clear  then,  the  present  problem  is  a 
cooperative  blend  of  the  two.  Just  as  conservation 
equations  for  the  mean  turbulent  stresses. 


can  be  derived  but  with  its  accompanying  closure 
problesis  (sec,  for  Instance,  (27)),  the  conservetlon 
equations  for  the  moduleted  turbulent  stresses, 
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can  also  be  similarly  obtained  (see,  for  instance, 
(28))  and  similar  problems  with  respect  to  closure 
also  follows.  Eddy  viscosity  models  for 


has  been  used  to  some  advantage  In  describing 
certain  features  of  the  large-scale  organized 
structure  (12),  (15),  (16),  (28),  however,  some 
progress  towards  removing  this  stigma  by  bringing 
In  the  active  participation  of  the  fine-grained 
turbulence  has  recently  been  made  In  a paper  which 
we  have  discussed  earlier  (24),  where  the  temporel 
mixing  region  problem  wee  discussed.  The  ideas 
follow  some  of  the  earlier  works  on  nonlinear  hydro- 
dynamic  stability  theory  for  laminar  shear  flows 
(26),  (29),  (30),  (31)  which  have  been  found  helpful 
in  explelnlng  the  large-acale  organized  structure 
In  developing  laminar  shear  flows.  In  this  paper, 
we  shall  present  our  work  on  the  spatially  develop- 
ing mixing  region  problem. 


ENERGY  INTEGRAL  DESCRIPTION 

The  considerations  alluded  to  in  our  dis- 
cussion of  Ref.  (24)  utilizes  techniques  of  non- 
linear hydrodynamic  stability  theory  while  bringing 
in  the  active  Interactions  with  the  fine-grained 
turbulence.  Such  a consideration  enables  one  to 
provide  a gross  description  of  the  Interaction 
processes  among  the  three  components  of  motion  and 
has  as  its  aim  the  derivation  of  the  respective 
coupled  "amplitude"  (or  "envelope")  equations, 
with  the  accompanying  shape  assumptions  for  the 
vertical  distributions  (24),  (26).  The  amplitude 
equations  are  derived  from  the  kinetic  energy 
equations  in  the  mean. 

For  the  spatially  developing  mixing  region, 
the  boundary  layer  approximations  are  applied  to 
the  mean  energy  equations  of  the  three  components 
of  flow.  The  application  of  such  approximations 
to  the  mean  flow  quantities  and  to  the  time- 
averaged  random  turbulence  quantities  hes  been 
the  usual  practice.  Such  approximations  to  the 
time-averaged  latge-scale  organized  structure 
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quantities  require  some  care.  In  the  previous 
worts  on  laminar  flows  such  approximations  were 
found  unnecessary  in  that  the  atreamvise  derivatives 
of  the  time-averaged  correlations  of  the  organized 
structure  could  be  evaluated  (29),  (31).  In  Ref. 
(31)  it  was  found  that  the  purely  dynamical  contri- 
butions were  unimportant  and  this  is  confirmed  in 
the  present  problem  also.  Rather  than  retaining 
these  terms,  they  are  neglected  at  the  outset  for 
simplicity.  The  appropriate  kinetic  energy 
equations  are  Chen 
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and  c is  the  rate  of  viscous  dlsslpstlon  of  the 
fine-grained  turbulence.  Viscous  diffusion  of  the 
kinetic  energies  snd  dlsslpstlon  of  the  mean  flow 
and  large-scale  structure  have  been  neglected.  In 
the  above,  che  respective  diffusion  and  production 
from  che  mean  flow  terms  are  self-explanatory  in 
che  two  fluctuation  energy  equations.  The  less 
familiar  average  energy  transfer  between  the  dis- 
parate scales  are  so  indicated. 

The  integral  kinetic  energy  equations  are  then 
r 
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with  the  upper  scream  V is  normalized  to  unity 
and  subsequently  the  lower  shear  stream  U taken 
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Implicit  in  the  above  is  that  the  mean  shear 
flow  is  two  dimensional,  similarly  with  the  large- 
scale  organized  structure  which,  in  the  subsequent 
shape  assumption,  will  be  taken  as  one  whose 
vorticity  axis  is  perpendicular  to  che  streamwlse 
direction.  The  recently  observed  spanwlac  period- 
icity (4)  is  not  accounted  here.  The  fine-grained 
turbulence,  of  course,  la  three  dimensional. 

The  closure  problem  is  obtained  in  the  dis- 
cussion of  the  shape  assumptions.  There,  some  of 
the  detailed  "microscopic"  calculations  of  the 
modulated  stresses  and  large-scale  structure,  prior 
to  averaging,  shall  be  done  approximately. 

Shape  Assumptions 

The  mean  flow  is  taken  to  be  of  the  form 
U - j (1  + tanh  n) 

where  n • y/6(x)  and  6(x)  is  the  x-dependenc 
shear  layer  thickness  which  will  characterize  the 
mean  flow,  to  be  determined  Jointly  with  the  large- 
scale  structure  and  the  fine-grained  turbulence 
problems. 

The  shape  functions  for  the  vertical  distri- 
bution of  the  time-averaged  fine-grained  turbulence 
quantities  are  guided  by  measurements  such  as  those 
of  Wygnanskl  and  Fiedler  (32).  The  slight  asymmetry 
of  the  profiles,  though  it  could  be  accomodated  here, 
is  not  taken  into  account.  The  Reynolds  shear 
stress  is  taken  to  be  of  the  form 
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and  is  related  to  the  turbulent  kinetic  energy  via 
the  constant  * 0.3.  thus  q • - u’v’/s,.  The 

factor  Jr  is  a normalizing  constant  so  that  fE(x) 
has  the  physical  interpretation  as  being  the  total 
kinetic  energy  of  the  fine-grained  turbulence  across 
a slice  of  the  entire  shear  layer.  E(x)  is  thus 
an  energy  density,  the  total  energy  per  unit  f(x) 
and  it  characterizes  the  turbulence.  Following 
Patankar  and  Spalding  (33),  and  Bradshaw.  Femss 
and  Atwell  (34),  the  dissipation  of 
the  fine-grained  turbulence  is  taken  to  be 

c - a;53/2/2a 

where  the  value  a,  - 1.5  will  be  used. 


11.4 


Following  earlier  works  (24),  (26),  (30), 
the  vertical  shape  of  the  large-scale  structure 
is  given  by  a toczl  linearized  thecrv  but  with 


ar.  amplitude  A(x)  which  is  to  be  Jointly  solved 
with  4(x)  anc  E(x), 


[U]  - A(x)e‘ifiT 

r q'(n;S) 

* 

-iap(n;£) 

)C.C. 

l J 

l J 

where  a is  the  dimensionless  complex  wave  number 
scaled  by  6(x),  6 is  the  dimensionless  physical 

frequency  scaled  also  according  to  the  local  shear 
layer  thickness  S(x)  and  is  real,  q is  the 
streamfunction  and  t’  refers  to  its  n-derlvative, 
c.c.  denotes  the  complex  conjugate.  Here  T is 
the  dimensionless  time  also  referred  to  6(x).  The 
local  shape  functions  are  coupled  to  and  diffused 
by  the  siodulated  turbulent  stresses.  The  structure 
of  these  stresses  warrants  their  shape  assumption 
tc  be  of  the  form  (24), 

~1£t 

- A(x)E(x)e  ri^(n;8)  c.c. 


Hence  the  vertical  shape  functions  q and  r^ 

forms  a coupled  local  linear  problem  and  in 
obtaining  their  equations  A is  taken  exponential 
in  x and  E constant  locally.  An  approximate 
method  of  solving  the  local  temporal  problem  is 
outlined  and  argued  in  (24)  and  can  be  carried  out 
similarly  for  the  spatial  problem  here.  Namely, 
the  "lnviscid"  form  of  £ is  first  solved  end 
then  used  in  the  algebraic  r^  equations  as 

known  functions.  Additional  closure  assumptions 
sre  required  for  use  in  the  equations  for  the 
modulated  stress  functions  r,,.  For  the  mean 

stresaes,  these  are  v - w and 

v '“*/(-  u’v')  * 2 from,  say,  (32).  For  the 
nodulated  turbulent  quantities  these  are 
patterned  after  closure  assumptions  for  the  mean 
quantities,  with  the  preasure  velocity  correlation 
taken  simply  as 


] 


• iff  -it  i 

T ij  rii  ij 

where  T is  the  time-scale  for  "return  to  iso- 
tropy" and  la  taken  as  T”1  • c8l)/3y  with 
c 4 1.445.  The  modulated  dissipation  is  taker 
to  be 

t S <t»  - C • j yi  $ 

where  8 * f,,/2.  Since  the  vertical  shape 
function  here1!*  to  be  given  by  the  local  linear 
theory,  the  form  for  i as  given  above  reflects 
a linearization  with  q there  to  be  replaced  by 
its  local  equilibrium  value  in  the  abaence  of 
the  large-scale  structure.  Then 


where  the  constant  t 4 0.24  (there  is,  of  course, 
no  confusion  with  the  same  symbol  used  for  the 
frequency) . This  is  similar  to  using  the  local 
homogeneous  shear  problem  as  a guide  (see.  for 
instance.  (24)).  These  additional  assumptions  are 
used  in  the  c,  r problem  which  enters  the  main 


overall  interaction  problem  for  d.  E and  A as 
quantities  occurring  under  "interaction"  integrals 
and  are  thus  not  considered  to  be  sensitive.  For 
completeness,  we  state  the  approximate  q,  r,, 

problem,  with  q governed  by  the  Rayleigh  equation 


(U  - c)(q"  - a2q)  - qU"  - 0 


with  appropriate  boundary  conditions,  and  r by 

* J \ 
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where  the  vertical  shape  functions  of  the  mean 

stresses  sre 


- 

- - 

R 

XX 

0.45 

R 

yy 

■ 

0.34 

R 

zz 

0.34 

-R 

*y 

_ 

0.17 

and  the  prime  on  the  R and  q denotes,  again, 

differentiation  with  respect  to  n.  Calculations 
for  q and  r^  are  ultimately  to  be  used  in  the 

integral  over  the  large-scale  structure  and 
turbulence  eneTgy  transfer  mechanism. 


VV»Y 


It  la  clear  that,  to  the  present  order  of  approxi- 
mations, such  a tern  is  proportional  to  ]a]2E  and 
the  "constants"  of  proportionality  are  functions  of 
4(x).  The  eigenfunctions  q are  normalized  so  as 
to  render  | A(x) | 2 to  be  the  corresponding  energy 
density  of  the  large-scale  organized  structure. 

Both  | A |2  and  E are  thus  similarly  interpreted. 


The  Interaction  Problem 


After  the  substitution  of  the  shape  functions 
into  the  energy  Integral  equations,  and  with  some 
straightforward  algebraic  dexterity,  we  obtain 


I.  ~ - I'  E - 1 (d ) - 1 A ; 2 

1 ox  rs  rs 


x3  2F " V ♦ - V3'2 


11.5 


12({)  ^dx1"  " ir.(S)'lAl2  - V«>-W*E 

with  Initial  condition*  6(0)  • 1,  E(0)  ■ E , 

o 

| A (0) | 2 • |a|2.  The  physical  interpretation  o i 

each  term  is  obvious  and  can  be  reverted  to  the 
original  conservation  equations. 

The  integrals  involving  only  the  mean  flow 
and  time-averaged  fine-grained  turbulence  quanti- 
ties are  constants, 

Ix  - j*  (1  + tanh  n)3dn 

+ i |”(1  + t»nh  n)jj  (1  * t*nh  n)2-  ljdn 

- - 0.202 


I„  • — i-  | (1  + t*nh  n)«  n dn  ■ 0.500 

3 lA  J- 


2 2 

I1  • — a*ch2n  e r dn  • 0.109 

r*  2^  J- 

' ft*  f. ' 


Th*  integral*  involving  eh#  l*rg#-»c*le  organized 
quantitle*  are  dependent  on  the  local  ahear  layer 
thlckneaa  6(x)  through  their  dependence  on  th* 
apectral  propertiea  via  local  frequency  parameter 
8(6).  The  large-acale  atructur*  production 
integral  la 

irg(6)  • iri(6)  “ j"  **ch2ndn 

where  Ei  denote*  th*  reel  part  and  t la  the 
complex  conjugate  of  ♦,  and  th*  large-acale 
structure  — fine-grained  turbulence  interaction 
integral  is 

Iwt(6)  ’ IwtU)  ' f Ri[rxx('1“*’)  + rxy(*"  + “2J) 

+ r (in*’)  dn 

Tf  J 

with  I >0.  Th#  detail*  of  th*  interpretation 

of  the  integrand  are  given  in  (24)  for  the  temporal 
problem,  and  will  not  be  repeated  here.  Finally, 


I2(S)  • lj(5)  * j jl  + j”  tanh  n-["]d 


+ |ar|b p dn 


and  is  very  nearly  conatant. 


The  spectral  properties  of  the  interaction 
integrals  associated  with  the  large-scale  structure 
are  shown  in  Figure  1.  They  are  numerically 
obtained  as  functions  of  6(x)  “ 6q6(x),  where  6o 

is  the  initial  cnaracteristic  (dimensionless) 
frequency  of  the  large-scale  structure,  B(0)  • co- 


Figure  1.  The  spectral  properties  of  the 
Interaction  integrals  involving  the  large-scale 
structure. 


The  advection  Integral  I,  is  very  nearly  constant 

and  it  will  be  approximated  by  0.5  during  the  sub- 
sequent calculations  since  ic  was  found  that  the 
problem  is  not  sensitive  to  the  variation  of  I,. 

The  Reynolds  stress  or  production  integral,  1 , 

peaks  at  about  £ * 0.21,  that  value  of  S which 
corresponds  to  the  maximal  amplification  in  the 
local  linear  theory.  The  large-acale  structure  — 
fine-grained  energy  cranafer  integral  I shown 

here  la  different  from  that  obtained  for  the 
temporal  problem  in  Ref.  (24)  where  the  spectral 
properties  were  plotted  as  functions  of  the  wave 
number.  There  the  local  turbulence  was  taken  to 
be  in  equilibrium  with  the  mean  flow  in  order  to 
isolate  and  illustrate  the  large-scale  structure 
and  fine-grained  turbulence  interaction  and  this 
rendered  the  modulated  stresses  to  become  singular 
as  the  local  linear  problem  approached  neutral. 
Here,  because  the  fine-grained  turbulence  is  not 
necessarily  in  equilibrium  with  the  mean  flow 
locally  no  such  singularity  occurred  as  6 -»  0.5. 
Calculations  of  the  properties  beyond  £ • 0.5 

were  not  obtained.  The  remaining  integral,  I , 

a 

is  shown  only  for  illustrative  purposes.  It  is 
associated  with  one  of  the  streamvise  large-scale 


< 

1 
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structure  diffusion  terns, 


and  the  enhanced  E (because  of  I A j 2 i becomes 
important,  the  growth  being  contribute..  bv  both  IX 
and  III. 


5(uf  + br  )/3x. 
xx  xy 

Since  I.  <<  I such  streamvise  diffusion  terms 
d vt 

contribute  very  little  to  the  main  features  of  the 
problem  and  they  are  neglected.  There  is  no  obvious 
rational  transformation  from  a temporal  to  a spatial 
problem  even  in  the  present  formulation. 

RESULTS 

With  the  large-scale  structure  interaction 
integrals  interpreted  as  functions  of  £ m Bo6(x) 

and  the  single  well  controlled  mode  la  specified 
to  be  £ , the  three  first  order  nonlinear  Inter- 
action equations  arc  Integrated  subjected  to  the 
initial  conditions  6(0)  • 1,  |a(0)|2  • | A | 2 . 

E(0)  ■ E . Energy  densities  corresponding  to 
-4 

about  10  Implies  root  mean  square  velocities  of 
about  a few  percent  of  the  free  stream  and  this 
seems  to  be  a reasonable  value  to  be  imposed 
upstream.  The  results  are  shown  as  functions  of 
x.  the  streamwlse  distance  normalized  by  the 
initial  mixing  layer  thickness . 

Physical  Mechanisms  for  the  Evolution  of  6 , 1 A I 2 , 
and  E 


In  Figure  2,  for  the  caae  of  6o  * 0.1,  two 

sets  of  curves  are  shown.  The  dashed  curve 
Indicates  the  shear  layer  growth  accomplished  by 
the  energy  transfer  from  the  mean  flov  to  the 
turbulence  energy  density  E alone.  Thus  the 
development  of  E is  due  to  the  imbalance  between 
production  from  the  mean  and  viscous  dissipation 
alone.  With  the  imposition  of  the  large-scale 
structure  the  mean  flow  spreading  becomes  enhanced 
as  is  the  turbulence  energy  density.  It  can  be 
seen  that  the  large-scale  structure  energy  density 
first  amplifies  and  subsequently  decays.  All  this 
is  quite  reminiscent  of  the  observations  of  Binder 
and  Favre-Marinet  (23).  In  order  to  obtain  physi- 
cal insight  into  the  behavior  shown  in  Figure  2, 
we  use  the  solution  in  a diagnostic  manner  by 
shoving  the  relevant  Interaction  terms  that  bring 
about  the  evolution  of  6,  E and  ] A | 2 . In  Figure  3 
we  illustrate  the  following  terms:  the  viscous 
dissipation  of  turbulence  energy 


shown  as  curve  I,  the  turbulence  production  from 
the  mean,  1^  E by  curve  II,  the  large-scale 

structure  production  from  the  mean  irj[A|2  by 

curve  III,  and  the  energy  transfer  between  the 
large-scale  and  turbulence  I I A 2E  bv  curve  TV. 

Initially  III  > IV  and  the  large-scale  structure 
amplifies  and  starts  tc  decay  when  III  < IV  after 
x 5 85.  Of  course.  III  itself  decays  because  of 
the  enhanced  spreading  caused  by  III  itself. 
Initially  when  ;A  - is  still  relatively  small, 
the  increase  of  E is  due  entirely  to  II  > 1. 

It  continues  tc  increase  and  becomes  much  enhanced 
aue  to  II  and  II  in  spite  of  the  fact  that  I > II 
after  x • 55.  The  spread  of  the  shear  layer  thick 
ness  6 is  significantly  enhanced  when  both  jA:' 


X 


Figure  2.  Streamvise  development  cf  6,  A 2 and 

E.  Solid  line:  £ -0.1.  A 2 - E - 10"; 

o o o 

Dotted  line:  £ * 0,  A **  • 0 and  E • 10"4. 

o o o 


X 

Figure  1.  Energv  transfer  mechanisms  ( • 0.1. 


The  Spectral  Dependence  of  the  Interaction 

Previously,  Che  interaction  was  shown  for 

o “0.1.  The  case  6 “ 0.21  is  one  for  which  the 

o o 

initial  large-scale  structure  is  most  amplified 
according  to  the  local  linear  theory  and  it  is  thus 
expected  that  the  large-scale  structure  at  this 
frequency  will  be  the  most  efficient  extractor  of 
energy  from  the  mean  flow  and  would  thus  attain 
the  highest  peak.  For  the  same  initial  conditions, 
the  development  of  | A | 2 for  6q  “ 0.1,  0.21  and 

0.3  are  shown  in  Figure  4.  As  is  expected,  the 
higher  frequency  large-scale  structures  peaked 
closer  upstream  while  lower  frequency  one*  further 

downstream  (12),  and  the  S “0.21  case  attained 

o 

the  highest  peak.  The  development  of  S,  | A | 2 and 
E is  shown  in  Figure  5 for  the  case  Bq  “ 0.21 

and  0.3,  and  these  are  to  be  compared  with  Figure  2 
for  the  same  Initial  conditions.  The  higher 
frequency  components  have  shorter  atreamwise  life- 
times; although  | A | 2 may  have  attained  higher 
peaks  than  the  low  frequency  modes  when  It  was 
most  amplified,  the  enhanced  E may  not  be  as 
vigorous  because  ‘ r depends  not  only  on  the  magni- 
tude of  j A j 2 so  on  the  lifetime  of  the 

large-scale  as  well. 


X 

Figure  4.  Spectral  effects  on  tne  streamvise 
cevelopment  of  large-scale  structure  energy 

aensity,  A - “ E - 10"".  (i  - 0.21  is  the 
o c o 

Initially  most  amplified  case.) 

controlling  the  Large-Scale  Structure 

In  addition  to  controlling  the  large-scale 
organized  structure  by  direct  imposition  of  definite 
spectral  modes  at  a giver,  strength,  the  use  of 
fine-grained  turbulence  to  control  its  development 


can  also  be  achieved  (12).  In  Figure  6 we  show 

the  8 “0.21  noae  with  fixed  A-  but  with 

o o 

different  Initial  conditions  for  Eq.  The  initially 

enhanced  fine-grained  turbulence  suppresses  the  sub- 
sequent development  of  the  large-scale  structure, 
curve  II,  and  in  so  doing,  it  achieves  a lower 
subsequent  enhancement  as  indicated  by  E... 

£0  * 0-21  £o*  0.3 


0 10  20  30  40  50  0 10  20  30  40 

x x 


Figure  5.  Spectral  effects  on  the  streamwise 
development  of  l,  A : and  E.  Solid  iine: 

i A ; 2 “ E ■ 10  6;  Dotted  line:  £ * 0.  ;a'*  • 0 

oo  o ‘ o 

and  Eq  “ 10  “.  (Compare  also  with  Figure  2.) 


Discussion 

The  spatial  nonlinear  interaction  problem 
has  been  obtained  by  consistently  interpreting  the 
interaction  integrals  involving  the  large-scale 
structure,  ifg  and  I , as  functions  of  i(x) 

through  their  dependence  on  B(x)  “ 8 6(x)  in  the 

o 

same  spirit  as  (29),  (30),  and  (31),  as  well  as 
(24).  In  the  temporal  problem  (24),  the  conse- 
quences of  freezing  these  integrals  as  constants 
were  further  explored.  The  problem  then  became 
analytically  tractable  so  as  to  enable  us  to 
obtain  equilibrium  states  in  the  | A ! 2 - E plane 
and  to  study  the  nature  of  the  singularity  about 
such  states.  It  was  found  that  the  |A|2  - E 
plane  describes  a stable  spiral  (24).  The  spiral 
behavior  was  sufficiently  pronounced  so  as  to 
resemble  "bursts"  or  renewed  large-scale  structure 
activity  after  repeated  cycles,  which  are  not 
actually  periodic,  of  amplification  and  decav. 

By  a similar  freezing  of  the  interaction  integrals 
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ir.  the  spatial  problem  here,  a spiral  behavior 
about  equilibrium  points  were  also  £ound  but  were 
not  as  pronounced  as  those  in  (24).  The  difference 
is  not  so  ouch  as  betveen  spatial  and  temporal 
problems,  but  that  in  (24)  the  local  turbulence 
behavior  was  taken  to  be  in  equilibrium  with  Che 

mean  and  this  demanded  a much  larger  I for 

wt 

those  modes  of  the  large-scale  structure  approach- 
ing the  neutral  condition.  In  the  present  problem, 
the  previous  simplifying  condition  of  local  turbu- 
lence equilibrium  with  the  mean  is  removed  and  this 

rendered  I to  be  finite  as  8 -*  0.5.  The 
wt 

results  for  the  frozen  interaction  integrals  are 
not  presented  here  because  of  space  limitations. 


numerical  simulation  could  also  be  attempted.  To 
this  end,  the  senior  author  of  the  present  paper 
recently  posed  the  following  problem. 

To  fix  our  ideas,  let  us  revert  tc  the 
temporal  mixing  layer  problem  of  horizontally 
nomogeneous  and  oppositely  directed  streams.  The 
large-scale  structure  is  horizontally  periodic  and 
develops  in  time.  The  "proper"  oynamical  dependent 
variable  is  the  sum.  of  the  mean  and  the  large-scale 
structure,  denoted  by  the  script  notation 

i'i(v.y.''-)  “ E(y.t)  + u1(x,y,t) 

P(x, y,t)  • P(y,t)  + p(x,y,t) 


X 


Figure  6.  The  effect  of  increasing  the  initial 
fine-grained  turbulence  energv  level  on  me 
development  of  large-scale  structure  and  turbulence 

energv  densities,  f - 0.21,  1:  A - * I - 10"; 

o , oo 

II:  A 1 - 10-',  E . io“ . 

o o 


THE  PROSPECTS  FOR  NUMERICAL  SIMULATION 


The  full  Kavler-Stokea  equations  is  then  split  into 
the  mean  and  organized  part  and  a turbulence  part 
according  to  + u^  and  p • P + p'.  Again 


denoting  the  conditional  average  by  <(.)>, 
have  for  the  dynamical  variablaa 
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end  the  coupled,  yet  uncloaed,  equation  for  the 
conditionally  averaged  streeaea 
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While  much  physical  intuition  could  be  gained 
about  the  interactions  between  large-scale  structure 
and  turbulence  in  a shear  flow  from  the  approximate 
energy  Integral  considerations,  much  detail  about 
the  local  structure  is  lacking.  The  analysis  of 
the  large-scale  structure  presented  here  resembles 
the  analysis  of  nonlinear  hydrodvnamic  stability 
theory  in  laminar  flow  (26),  (29),  (30),  (31). 
Complementary  to  such  analysis  is  the  numerical 
simulation  of  f inite- amplitude  instabilities  in 
laminar  flow  using  the  full  nonsteadv  Navier- 
Stokes  equations  (see,  for  instance,  (35),  (36)). 

For  the  large-scale  organized  structure  ir.  a 
turbulent  shear  flow,  a similar  but  more  difficult 


dissipation 


The  viscous  diffusion  effects  will  subsequentlv  be 
neglected.  A possible  closure  scheme  is  obtainet 
bv  a direct  extension  of  Launder  et  al.'s  model 
(27)  for  the  mean  stresses  to  the  present  condi- 
tionally averaged  stresses  »u!u!>.  Again, 

reverting  to  the  temporal  mixing  layer  problem 
with  a two-dimensional  large-scale  structure, 
such  a model  then  provides  eight  equations  for 
the  eight  unknowns  l,  V,  F and  *u'v'»,  <u'J>, 
<v'->,  and  <c>.  (Or,  following  Amsdet  and 

Harlov  (35),  the  three  dynamical  variables  are 


replaced  by  Che  vortlcity  and  streamfunction. ) 

!'  Is  horizontally  periodic  as  are  the  condition- 
ally averaged  stresses  <u^uj>  and  these  decay  far 

from  the  shear  layer.  Originally  there  existed  the 
scan  mixing  layer  with  a set  of  consistent  and  com- 
patible mean  stresses  and  at  t ■ tQ,  a large-scale 

organized  structure  is  imposed  so  chat 


^(x.y.to) 


= U (y,t  ) + a (x.y.t  ), 

10  1 0 
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1357. 

5 Wlllmarth,  K.  W.,  "Structure  of  Turbulence 
in  Boundary  Layers,"  Aovances  in  Applied  Mechanics. 
Vol.  15,  Academic  Press,  New  York,  1975,  pp.  159- 
254. 


profile  and  t^Cx.y,^)  taken  as  zero. 

Ut  * h + ^i 


«u[u*»(x,y,to)  5 «iu'(y.t0)  + f^fx.y,^). 

Hence,  the  "initialization"  calls  for  a compatible 
coaxlstance  of  the  mean  problem  l^fy,^)  and 

uTuT(y,t  ).  The  impoaad  0, (x,y,t  ) is  then  made 

1 J O 10 

compatible  with  Ui(y,tQ),  being,  say,  the  solution 
to  the  Rayleigh  equation  with  U(y,tn)  as  the  mean 

The  problem 

yields  the  total  quantities,  U,  2 0,  + u,  and 
<u^uj>  i u^u^  + fy,  with  the  mean  obtainable  by 

horizontal  averaging  and  the  large-scale  structure 
quantities  by  a subtraction  of  the  mean  from  the 
total  quantities.  Of  great  interest,  for  instance, 
is  the  use  of  the  numerical  solution  for  diagnosing 
the  energy  transfer  between  the  large-scale  struc- 
ture and  che  fine-grained  turbulence, 

among  others.  The  numerical  problem  here  is 
certainly  compounded  compared  to  the  relatively 
more  straightforward  laminar  problem.  Our 
colleague,  T.  B.  Gatski,  is  courageously  working 
on  this  problem  and  the  results,  when  available, 
would  be  of  great  interest  towards  physically 
diagnosing  the  interaction  problem.  The  extension 
of  the  temporally  developing  to  the  spatially 
developing  mixing  layer  problem  is  certainly 
possible.  However,  much  of  tha  physics  about  the 
Interaction  could  be  gained  from  the  simpler, 
temporal  problem  first  before  attempting  che  more 
difficult  spatial  problem. 
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